library(tidyverse)
library(tuneR)
library(devtools)
library(ggplot2)
library(tsfeatures)
library(class)
library(cvTools)
library(randomForest)
library(GEOquery) 
library(R.utils)
library(reshape2)
library(limma)
library(dplyr)
library(e1071)
library(DT)
library(viridis)
library(plotly)
library(scales)
gse36059 = getGEO("GSE36059")[[1]]
Found 1 file(s)
GSE36059_series_matrix.txt.gz
Using locally cached version: C:\Users\lmcca\AppData\Local\Temp\RtmpC6OZAY/GSE36059_series_matrix.txt.gz
Warning in for (getMethodName in getMethodNames) { :
  closing unused connection 5 (https://ftp.ncbi.nlm.nih.gov/geo/series/GSE36nnn/GSE36059/matrix/)
Warning in .Internal(lapply(X, FUN)) :
  closing unused connection 4 (https://ftp.ncbi.nlm.nih.gov/geo/series/GSE36nnn/GSE36059/matrix/)
Warning in .Internal(lapply(X, FUN)) :
  closing unused connection 3 (https://ftp.ncbi.nlm.nih.gov/geo/series/GSE48nnn/GSE48581/matrix/)
Using locally cached version of GPL570 found here:
C:\Users\lmcca\AppData\Local\Temp\RtmpC6OZAY/GPL570.soft.gz 
gse48581 = getGEO("GSE48581")[[1]]
Found 1 file(s)
GSE48581_series_matrix.txt.gz
Using locally cached version of GPL570 found here:
C:\Users\lmcca\AppData\Local\Temp\RtmpC6OZAY/GPL570.soft.gz 
gse129166 = getGEO("GSE129166")[[1]]
Found 1 file(s)
GSE129166_series_matrix.txt.gz
Using locally cached version of GPL570 found here:
C:\Users\lmcca\AppData\Local\Temp\RtmpC6OZAY/GPL570.soft.gz 
gse36074_f = fData(gse36059)
gse36074_f
gse36074_f = fData(gse48581)
gse36074_f
gse129166_f = fData(gse129166)
gse129166_f
gse36074_p = pData(gse36059)
gse36074_p
gse48581_p = pData(gse48581)
gse48581_p
gse129166_p = pData(gse129166)
gse129166_p
gse36074_p$diagnosis = ifelse(gse36074_p$characteristics_ch1 == "diagnosis: non-rejecting", 0, 1)
gse48581_p$diagnosis = ifelse(gse48581_p$characteristics_ch1.1 == "diagnosis (tcmr, abmr, mixed, non-rejecting, nephrectomy): non-rejecting", 0, 1)
gse129166_p$diagnosis = ifelse((gse129166_p$characteristics_ch1.1 == "tcmr (no: 0_borderline:1_TCMR:2): 0") & (gse129166_p$characteristics_ch1.2 == "abmr (no: 0_Yes:1): 0"), 0, 1)
gse36074_p$exact_diagnosis = ifelse(gse36074_p$characteristics_ch1 == "diagnosis: non-rejecting", 0, ifelse(gse36074_p$characteristics_ch1 =="diagnosis: ABMR", 1, ifelse(gse36074_p$characteristics_ch1 == "diagnosis: TCMR", 2, 3)))

gse48581_p$exact_diagnosis = ifelse(gse48581_p$characteristics_ch1 == "diagnosis (tcmr, abmr, mixed, non-rejecting, nephrectomy): non-rejecting", 0, ifelse(gse48581_p$characteristics_ch1 =="diagnosis (tcmr, abmr, mixed, non-rejecting, nephrectomy): ABMR", 1, ifelse(gse48581_p$characteristics_ch1 == "diagnosis (tcmr, abmr, mixed, non-rejecting, nephrectomy): TCMR", 2, 3)))

gse129166_p$exact_diagnosis = ifelse((gse129166_p$characteristics_ch1.1 == "tcmr (no: 0_borderline:1_TCMR:2): 0") & (gse129166_p$characteristics_ch1.2 == "abmr (no: 0_Yes:1): 0"), 0, ifelse((gse129166_p$characteristics_ch1.1 != "tcmr (no: 0_borderline:1_TCMR:2): 0") & (gse129166_p$characteristics_ch1.2 != "abmr (no: 0_Yes:1): 0"), 3, ifelse((gse129166_p$characteristics_ch1.1 != "tcmr (no: 0_borderline:1_TCMR:2): 0") & (gse129166_p$characteristics_ch1.2 == "abmr (no: 0_Yes:1): 0"), 2, 1)))
gse36074_p
gse48581_p
gse129166_p

Main outcomes are non-rejecting, TCMR (acute T-cell–mediated rejection), ABMR (anti-donor antibody-mediated rejection), and MIXED. Also has nephrectomy.

gse36074_ex = data.frame(t(exprs(gse36059)))
gse36074_ex
gse48581_ex = data.frame(t(exprs(gse48581)))
gse48581_ex
gse129166_ex = data.frame(t(exprs(gse129166)))
gse129166_ex
start = 1
stop = 100
boxplot(gse36074_ex[start:stop])

boxplot(gse48581_ex[start:stop])

boxplot(gse129166_ex[start:stop])

LS0tDQp0aXRsZTogIkxpYW0tUmVzZWFyY2gtMi1JbnRlZ3JhdGlvbiINCm91dHB1dDogaHRtbF9ub3RlYm9vaw0KLS0tDQoNCmBgYHtyfQ0KbGlicmFyeSh0aWR5dmVyc2UpDQpsaWJyYXJ5KHR1bmVSKQ0KbGlicmFyeShkZXZ0b29scykNCmxpYnJhcnkoZ2dwbG90MikNCmxpYnJhcnkodHNmZWF0dXJlcykNCmxpYnJhcnkoY2xhc3MpDQpsaWJyYXJ5KGN2VG9vbHMpDQpsaWJyYXJ5KHJhbmRvbUZvcmVzdCkNCmxpYnJhcnkoR0VPcXVlcnkpIA0KbGlicmFyeShSLnV0aWxzKQ0KbGlicmFyeShyZXNoYXBlMikNCmxpYnJhcnkobGltbWEpDQpsaWJyYXJ5KGRwbHlyKQ0KbGlicmFyeShlMTA3MSkNCmxpYnJhcnkoRFQpDQpsaWJyYXJ5KHZpcmlkaXMpDQpsaWJyYXJ5KHBsb3RseSkNCmxpYnJhcnkoc2NhbGVzKQ0KYGBgDQoNCmBgYHtyfQ0KZ3NlMzYwNTkgPSBnZXRHRU8oIkdTRTM2MDU5IilbWzFdXQ0KZ3NlNDg1ODEgPSBnZXRHRU8oIkdTRTQ4NTgxIilbWzFdXQ0KZ3NlMTI5MTY2ID0gZ2V0R0VPKCJHU0UxMjkxNjYiKVtbMV1dDQpgYGANCg0KDQoNCmBgYHtyfQ0KZ3NlMzYwNzRfZiA9IGZEYXRhKGdzZTM2MDU5KQ0KZ3NlMzYwNzRfZg0KZ3NlMzYwNzRfZiA9IGZEYXRhKGdzZTQ4NTgxKQ0KZ3NlMzYwNzRfZg0KZ3NlMTI5MTY2X2YgPSBmRGF0YShnc2UxMjkxNjYpDQpnc2UxMjkxNjZfZg0KYGBgDQoNCmBgYHtyfQ0KZ3NlMzYwNzRfcCA9IHBEYXRhKGdzZTM2MDU5KQ0KZ3NlMzYwNzRfcA0KZ3NlNDg1ODFfcCA9IHBEYXRhKGdzZTQ4NTgxKQ0KZ3NlNDg1ODFfcA0KZ3NlMTI5MTY2X3AgPSBwRGF0YShnc2UxMjkxNjYpDQpnc2UxMjkxNjZfcA0KYGBgDQoNCg0KYGBge3J9DQpnc2UzNjA3NF9wJGRpYWdub3NpcyA9IGlmZWxzZShnc2UzNjA3NF9wJGNoYXJhY3RlcmlzdGljc19jaDEgPT0gImRpYWdub3Npczogbm9uLXJlamVjdGluZyIsIDAsIDEpDQpnc2U0ODU4MV9wJGRpYWdub3NpcyA9IGlmZWxzZShnc2U0ODU4MV9wJGNoYXJhY3RlcmlzdGljc19jaDEuMSA9PSAiZGlhZ25vc2lzICh0Y21yLCBhYm1yLCBtaXhlZCwgbm9uLXJlamVjdGluZywgbmVwaHJlY3RvbXkpOiBub24tcmVqZWN0aW5nIiwgMCwgMSkNCmdzZTEyOTE2Nl9wJGRpYWdub3NpcyA9IGlmZWxzZSgoZ3NlMTI5MTY2X3AkY2hhcmFjdGVyaXN0aWNzX2NoMS4xID09ICJ0Y21yIChubzogMF9ib3JkZXJsaW5lOjFfVENNUjoyKTogMCIpICYgKGdzZTEyOTE2Nl9wJGNoYXJhY3RlcmlzdGljc19jaDEuMiA9PSAiYWJtciAobm86IDBfWWVzOjEpOiAwIiksIDAsIDEpDQpgYGANCg0KYGBge3J9DQpnc2UzNjA3NF9wJGV4YWN0X2RpYWdub3NpcyA9IGlmZWxzZShnc2UzNjA3NF9wJGNoYXJhY3RlcmlzdGljc19jaDEgPT0gImRpYWdub3Npczogbm9uLXJlamVjdGluZyIsIDAsIGlmZWxzZShnc2UzNjA3NF9wJGNoYXJhY3RlcmlzdGljc19jaDEgPT0iZGlhZ25vc2lzOiBBQk1SIiwgMSwgaWZlbHNlKGdzZTM2MDc0X3AkY2hhcmFjdGVyaXN0aWNzX2NoMSA9PSAiZGlhZ25vc2lzOiBUQ01SIiwgMiwgMykpKQ0KDQpnc2U0ODU4MV9wJGV4YWN0X2RpYWdub3NpcyA9IGlmZWxzZShnc2U0ODU4MV9wJGNoYXJhY3RlcmlzdGljc19jaDEgPT0gImRpYWdub3NpcyAodGNtciwgYWJtciwgbWl4ZWQsIG5vbi1yZWplY3RpbmcsIG5lcGhyZWN0b215KTogbm9uLXJlamVjdGluZyIsIDAsIGlmZWxzZShnc2U0ODU4MV9wJGNoYXJhY3RlcmlzdGljc19jaDEgPT0iZGlhZ25vc2lzICh0Y21yLCBhYm1yLCBtaXhlZCwgbm9uLXJlamVjdGluZywgbmVwaHJlY3RvbXkpOiBBQk1SIiwgMSwgaWZlbHNlKGdzZTQ4NTgxX3AkY2hhcmFjdGVyaXN0aWNzX2NoMSA9PSAiZGlhZ25vc2lzICh0Y21yLCBhYm1yLCBtaXhlZCwgbm9uLXJlamVjdGluZywgbmVwaHJlY3RvbXkpOiBUQ01SIiwgMiwgMykpKQ0KDQpnc2UxMjkxNjZfcCRleGFjdF9kaWFnbm9zaXMgPSBpZmVsc2UoKGdzZTEyOTE2Nl9wJGNoYXJhY3RlcmlzdGljc19jaDEuMSA9PSAidGNtciAobm86IDBfYm9yZGVybGluZToxX1RDTVI6Mik6IDAiKSAmIChnc2UxMjkxNjZfcCRjaGFyYWN0ZXJpc3RpY3NfY2gxLjIgPT0gImFibXIgKG5vOiAwX1llczoxKTogMCIpLCAwLCBpZmVsc2UoKGdzZTEyOTE2Nl9wJGNoYXJhY3RlcmlzdGljc19jaDEuMSAhPSAidGNtciAobm86IDBfYm9yZGVybGluZToxX1RDTVI6Mik6IDAiKSAmIChnc2UxMjkxNjZfcCRjaGFyYWN0ZXJpc3RpY3NfY2gxLjIgIT0gImFibXIgKG5vOiAwX1llczoxKTogMCIpLCAzLCBpZmVsc2UoKGdzZTEyOTE2Nl9wJGNoYXJhY3RlcmlzdGljc19jaDEuMSAhPSAidGNtciAobm86IDBfYm9yZGVybGluZToxX1RDTVI6Mik6IDAiKSAmIChnc2UxMjkxNjZfcCRjaGFyYWN0ZXJpc3RpY3NfY2gxLjIgPT0gImFibXIgKG5vOiAwX1llczoxKTogMCIpLCAyLCAxKSkpDQoNCmBgYA0KDQoNCmBgYHtyfQ0KZ3NlMzYwNzRfcA0KZ3NlNDg1ODFfcA0KZ3NlMTI5MTY2X3ANCmBgYA0KTWFpbiBvdXRjb21lcyBhcmUgbm9uLXJlamVjdGluZywgVENNUiAoYWN1dGUgVC1jZWxs4oCTbWVkaWF0ZWQgcmVqZWN0aW9uKSwgQUJNUiAoYW50aS1kb25vciBhbnRpYm9keS1tZWRpYXRlZCByZWplY3Rpb24pLCBhbmQgTUlYRUQuIEFsc28gaGFzIG5lcGhyZWN0b215Lg0KDQoNCmBgYHtyfQ0KZ3NlMzYwNzRfZXggPSBkYXRhLmZyYW1lKHQoZXhwcnMoZ3NlMzYwNTkpKSkNCmdzZTM2MDc0X2V4DQpnc2U0ODU4MV9leCA9IGRhdGEuZnJhbWUodChleHBycyhnc2U0ODU4MSkpKQ0KZ3NlNDg1ODFfZXgNCmdzZTEyOTE2Nl9leCA9IGRhdGEuZnJhbWUodChleHBycyhnc2UxMjkxNjYpKSkNCmdzZTEyOTE2Nl9leA0KYGBgDQpgYGB7cn0NCnN0YXJ0ID0gMQ0Kc3RvcCA9IDEwMA0KYm94cGxvdChnc2UzNjA3NF9leFtzdGFydDpzdG9wXSkNCmJveHBsb3QoZ3NlNDg1ODFfZXhbc3RhcnQ6c3RvcF0pDQpib3hwbG90KGdzZTEyOTE2Nl9leFtzdGFydDpzdG9wXSkNCmBgYA0KDQoNCg0KDQpgYGB7cn0NCm1lcmdlKGdzZTM2MDc0X2V4LCBnc2U0ODU4MV9leCwgYWxsID0gVFJVRSkNCmBgYA0KDQoNCg0KDQoNCg0KDQoNCg0KDQoNCg0KDQoNCg0KDQoNCg0KDQoNCg0KDQoNCg0KDQo=